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ABSTRACT 

We present a statistical study of a large, homogeneously analyzed sample of narrow- 
line Seyfert 1 (NLS1) galaxies, accompanied by a comparison sample of broad-line 
Seyfert 1 (BLS1) galaxies. Optical emission-line and continuum properties are sub- 
jected to correlation analyses, in order to identify the main drivers of active galactic 
nuclei (AGN) correlation space, and of NLS1 galaxies in particular. For the first time, 
we have established the density of the narrow-line region as a key parameter in Eigen- 
vector 1 space, as important as the Eddington ratio L/L^d- This is important because 
it links the properties of the central engine with the properties of the host galaxy; i.e., 
the interstellar medium (ISM). We also confirm previously found correlations involv- 
ing the line width of H/3, and the strength of the Fe II and [O III] A5007 emission 
lines, and we confirm the important role played by L/Lgdd m driving the properties of 
NLS1 galaxies. A spatial correlation analysis shows that large-scale environments of the 
BLS1 and NLS1 galaxies of our sample are similar. If mergers are rare in our sample, 
accretion-driven winds on the one hand, or bar-driven inflows on the other hand, may 
account for the strong dependence of Eigenvector 1 on ISM density. 

Subject headings: galaxies: ISM - galaxies: active - galaxies: Seyfert - quasars: emis- 
sion lines 

1. Introduction 

Narrow-line Seyfert 1 (NLS1) galaxies as active galactic nuclei (AGN) with the narrowest 
Balmer lines from the broad-line region (BLR) and the strongest Fe II emission, cluster at one 
extreme end of AGN correlation space. It is expected that such correlations provide some of the 
strongest constraints on, and new insights in, the physical conditions in the centers of AGN and 
the prime drivers of activity, and the study of NLS1 galaxies is therefore of particular interest (see 
Komossa 2008 for a review). 

A number of key results have arisen from principal component analyses (PC A). Applied to 
NLS1 and broad-line Seyfert 1 (BLS1) galaxies, the strongest correlations are among optical AGN 
properties, and involve the width of the H/3 line and the strength of the [O III] line and Fe II 
complex. Others include the steepness of the X-ray spectrum and the asymmetry of the CIV 
emission line (e.g. Boroson Sz Green 1992, Laor et al. 1997, Brandt 1999, Marziani et al. 2001, 
Sulentic et al. 2000, 2002, Boroson 2002, Grupe 2004, Sulentic et al. 2007, 2008). 



These multi-wavelength observations (and their interpretation or modeling) hint at smaller 
black hole masses in NLS1 galaxies, and as such their black holes represent an important link with 
the elusive intermediate mass black holes, which have been little studied so far. Accreting likely 
at very close to the maximum allowed values, NLS1 galaxies are important test-beds of accretion 
models. 

Given their low black hole masses, mergers of NLS1 galaxies are prime sources of gravitational 
wave emission within the LISA sensitivity band (e.g., Centrella 2010). 

Independent samples are of importance when assessing the robustness of correlation analyses, 
and to increase correlation space (Grupe et al. 1999, Veron-Cetty et al. 2001, Williams et al. 
2002, Zhou et al. 2006, Grupe et al. 2010). We have studied a sample of 94 AGN, including 55 
NLS1 galaxies and 39 BLS1 galaxies. Ours is a large, homogeneously analyzed sample, and we 
add new emission-line measurements (particularly, the density-sensitive [S II] ratio) to correlation 
analyses. The comparison sample of BLS1 galaxies enables us to study trends across the whole 
AGN population; and to distinguish between physical processes operating only within the NLS1 
population, and those which are relevant for AGN as a class. 

This paper is the fourth in a sequence. The first one concentrated on the differences in the 
narrow-line region (NLR) density of BLS1 and NLS1 galaxies (Xu et al. 2007, hereafter X07), 
the second demonstrated that NLS1 and BLS1 galaxies follow the same Mbh — & relation if "blue 
outliers" are removed from the sample (Komossa k, Xu 2007, hereafter KX07), and the third 
discussed the properties of Seyfert galaxies with extreme outflows ("blue outliers"; Komossa et 
al. 2008, hereafter K08). In this fourth paper of the sequence, we focus on sample properties 
and correlation analyses. This paper is organized as follows. In Section 2 we describe the sample 
selection and the methods of data analysis. Section 3 reports the results obtained from the emission- 
line and continuum measurements. In Section 4 and 5 results from our correlation analyses and 
cross-correlation analysis are given, which are then discussed and interpreted in Section 6. We 
conclude with a summary in Section 7. 

We use the terms NLS1 galaxies and BLS1 galaxies collectively for high-luminosity and low- 
luminosity objects, i.e., jointly for Seyfert galaxies and quasars. Throughout this paper, a cosmology 
with Hq = 70 km s _1 Mpc -1 , = 0.3 and Q\ = 0.7 is adopted. 

2. Data analysis 
2.1. The sample 

Our sample consists of NLS1 galaxies from the catalog of Veron-Cetty & Veron (2003), to 
which we added a comparison sample of BLS1 galaxies from Boroson (2003) at z < 0.3, first 
presented by X07. All galaxies have been observed in the course of the Sloan Digital Sky Survey 
(SDSS) (Data Release 3; Abazajian et al. 2005) and have detectable low-ionization emission lines 



(in particular, [S II] AA6716, 6731 is always present with S/N > 5). The sample of NLS1 galaxies 
has similar redshift and absolute magnitude distribution as the BLS1 galaxy sample. The original 
sample selection, data preparation, and data analysis methods are described in detail by X07. 
Here, we briefly summarize the procedures in Sect. 2.2 and 2.3. For the first time, we present our 
measurements for the full sample in table form (Table 1). While the focus of our study is on NLS1 
galaxies, having a comparison sample of BLS1 galaxies is important. Our BLS1 sample contains 
approximately equal numbers as the NLS1 sample. While this is the best approach for comparison 
purposes, we also note that in observations, NLS1 galaxies make only approximately 20% of the 
Seyfert 1 population. 



2.2. Decomposition of starlight and nuclear continuum 

All SDSS spectra were corrected for Galactic extinction, the continuum was decomposed into 
host galaxy and AGN components (Lu et al. 2006), and then we subtracted the starlight component, 
nuclear continuum and the Fe II complexes from the spectra. More precisely, the spectra were 
decomposed into the following four components (see X07 for details): (1) A starlight component 
modeled by 6 synthesized galaxy templates, which were built from the synthetic spectral library of 
Bruzual & Chariot (2003). The details of the algorithm were described in Lu et al. (2006). (2) A 
power-law continuum to represent the AGN continuum. (3) A Balmer continuum generated in the 
same way as Dietrich et al. (2003). (4) An Fe II template given by Veron-Cetty & Veron (2004). 
The final multicomponent fit including the Fe II complexes was then subtracted from the observed 
spectrum and the emission line properties were measured in the way described in Sect. 2.3. 



2.3. Emission line fits 

Emission line profiles of the galaxies were fit with Gaussians and/or Lorentzians using the IRAF 
package SPECFIT (Kriss 1994). Measured FWHMs were corrected for instrumental broadening. 
Complex emission lines, which cannot be represented by a single Gaussian profile, were decomposed 
into multiple components. The Balmer lines were decomposed into a narrow and a broad component 
representing emission from the NLR and BLR, respectively. The narrow core was fit employing a 
single Gaussian profile with FWHM fixed to that determined for [S II] AA6716, 6731. As in previous 
work (e.g., Rodriguez-Ardila et al. 2000, Dietrich et al. 2005, X07, Mullaney & Ward 2008), the 
broad component was fit by using a combination of two Gaussian profiles, or alternatively a single 
Lorentzian profile (e.g., Veron-Cetty et al. 2001, Sulentic et al. 2002). For the approach employing 
the two-Gaussian profile fit, the final width of the broad-line emission, H/?b, is determined as the 
FWHM of the sum of the two Gaussians. We do not assign a physical meaning to the two separate 
broad components; they merely serve as a mathematical description of the complex shapes of the 
broad lines. While comparable and equally reasonable results can be achieved with both single 
Lorentzian and two Gaussian profiles for most broad components of the Balmer lines of NLS1 



galaxies, generally no acceptable fit is possible when employing the Lorentzian profiles to fit BLS1 
galaxies. Among the NLS1 galaxies, line width measurements between Lorentzian and Gaussian fits 
typically agree within 20%, Lorentzian profiles always resulting in smaller FWHMs, as expected. 
We finally used the results of the multi- Gaussian fits (Table 1) for the NLS1 classification and also 
for further correlation analysis. This allows a direct comparison with the BLS1 control sample, and 
also with several previous studies (e.g. Boroson & Green 1992, Grupe et al. 1999, Vaughan et al. 
2001). 

Apart from the Balmer lines, [O III] shows a complex profile. The total [O III] profile, [O IIl] to tai 
was decomposed into two Gaussian components: a narrow core and a broad base. The narrow core 
of [O III] is referred to as [O III] C . Measurements of the FWHM and blueshift of [O III] (Table 1) 
refer to the core of [O III], unless noted otherwise. The velocity shift of [O III] was measured relative 
to [S II]. We use positive velocity values to refer to blueshifts, negative ones for redshifts. All other 
forbidden lines were well fit with single Gaussian profiles. 

The strength of the Fe II emission, R4570, was measured as the ratio of the flux of the Fe II 
complex between the rest wavelengths 4434 and 4684 A to that of total H/3 emission. The strength 
of the [O III] emission, R5007, was measured as the flux ratio of total [O III] A5007 over total H/3 
emission. "Total" H/3 emission refers to the sum of broad and narrow component. 

2.4. A note on NLS1 classification 

After re-classification based on spectral emission- line fitting, we have 39 BLS1 and 55 NLS1 
galaxies in our sample. The standard classification criterion of NLS1 galaxies is according to 
FWHM(H/3 b ) < 2000 km s" 1 (Goodrich 1989) which is most commonly applied. Such NLS1 
galaxies usually come with weak [O IIl]A5007/H/3 and strong Fe II/H/3 which form part of the 
classification criteria (Osterbrock & Pogge 1985, Goodrich 1989, Veron-Cetty et al. 2001). We 
note in passing that Sulentic et al. (2001) rather suggest to use FWHM(H/3b) = 4000 km s _1 
as 'dividing' value, because they find AGN properties to change more significantly at that value. 
Furthermore, a few authors pointed out that NLS1 classification should incorporate a luminosity 
dependence of the FWHM(H/3 b ) cut-off (e.g., Laor 2000, Veron-Cetty et al. 2001, Shemmer et al. 
2004). Netzer & Trakhtenbrot (2007) suggest to base NLS1 classification on the Eddington ratio, 
and use L/LEdd > 0.25 for NLS1 galaxies. In the current study, we continue to use the classical 
FWHM cut-off of FWHM(H/3 b ) < 2000 km s" 1 which is still most commonly used. Results (Table 1) 
are reported in dependence of FWHM(H/3b) and L/L^d, so that other cut-off values could be easily 
applied to our sample. 



3. Emission-line and AGN properties 



Results from our emission- line and continuum measurements are listed in Table 1. Relevant 
parameters of the continua and emission lines were used to derive some AGN parameters, including 
black hole masses and Eddington ratios (Figure 1). We describe below how we determined these 
parameters. 

Assuming that the motion of the BLR clouds is virialized (e.g., Wandel et al. 1999), the 
black hole mass can be estimated as A/bh = Rbirv 2 / G . The velocity v of the BLR clouds can be 
estimated from the FWHM as v = V3/2 x FWHM, by assuming an isotropic cloud distribution. 
There is only few reverberation mapping results for NLS1 galaxies so far (Peterson et. al. 2000, 
Bentz et al. 2009). We assumed that the relation between the radius of the BLR and the optical 
luminosity, established from reverberation mapping for nearby Seyfert galaxies (e.g., Peterson et 
al. 2004; Kaspi et al. 2005) is indeed applicable to NLS1 galaxies (e.g., Komossa et al. 2006). 

We estimated the black hole masses of our NLS1 and BLS1 galaxies using the Rblr~^L^iqq 
relation given by Kaspi et al. (2005). The luminosity L5100 was derived from the SDSS photometr^. 
The PSF g and r magnitudes (Galactic extinction corrected) were used to determine the continuum 
slope (e.g., Wu & Liu 2004) and the flux density at 5100 A rest frame. The FWHM of the broad 
component from the multi- Gaussian fit of H/3, FWHM(H/3b) (Table 1), was used for the black hole 
mass estimates reported in Table 1. An error in black hole mass is typically 0.5 dex, which arises 
from the use of single-epoch data (e.g., Vestergaard 2004), and uncertainties in H/3 decomposition 
and in host galaxy contribution. 

The NLS1 galaxies of our sample have, on average, smaller black hole masses than BLS1 
galaxies. The estimated black hole masses of NLS1 galaxies range from log -Mbh,nlsi = 5.7 M to 
7.3 M Q with an average value of 6.5 M Q , while BLS1 galaxies range from log -Mbh.blsi = 6.5 M to 
8.4 Mq with an average value of 7.2 Mq (Figure 1). 

Eddington ratios, represented by L/L^dd, were calculated from the deduced black hole masses, 
according to ^Edd = 1-3 10 38 Mbr/Mq erg s -1 . For bolometric luminosity correction, we assumed 
-^boi = 9 AL5100 (Kaspi et al. 2000). As shown previously, NLS1 galaxies exhibit, on average, higher 
Eddington ratios than BLS1 galaxies. Eddington ratios of the NLS1 galaxies of our sample range 
from log £/£Edd,NLSi = ~~ 0.6 to 0.3 with an average value of —0.1, while BLS1 galaxies range from 



The SDSS photometry calibration is more robust given that the imaging data is taken on moonless photometric 
nights under good seeing conditions (Hogg et al. 2001), while spectroscopy is done on those nights that are not 
photometric, or with seeing worse than 1.7" FWHM, or with a moderate amount of moon (e.g., Abazajian et al. 
2003). Moreover, the SDSS spectrophotometric fluxes before Data Release 6 were tied to fiber magnitudes that led 
to an under-estimation in the spectrophotometric scale (Adelman-McCarthy et al. 2008) . Magnitudes used to derive 
black hole masses include a contribution from host galaxies. However, most of our objects are AGN dominated and 
thus the contamination from host galaxies is small. We compared our fluxes derived from the SDSS photometry to 
the spectroscopic fluxes and found that differences of flux scale between the two calibrations can reach at maximum 
35%. 



log i/^Edd,BLSi = —1-4 to —0.3 with an average value of —0.8. 

The distributions of black hole masses and Eddington ratios are shown in Figure la and lb, 
respectively. We have also applied a Kolmogorov-Smirnov (K-S) test, in order to confirm that 
the two distributions for NLS1 and BLS1 galaxies are significantly different. Indeed, for both the 
distributions, we find a K-S probability of P < 10~ 4 that the two samples are drawn from the same 
parent population. For comparison, we also plot histograms of other properties (Figure 2); some of 
these reflecting trends we reported previously. (1) The distribution of FWHM([S II]) (Figure 2a) 
shows narrower [S II] in NLS1 galaxies than BLS1 galaxies (KX07). (2) The distribution of [S II] 
intensity ratios and the NLR density (Figure 2b and 2c) demonstrate our discovery of a zone of 
avoidance in density, in the sense that BLS1 galaxies avoid low densities, while NLS1 galaxies show 
a larger scatter in density, including a significant number of objects with low densities (X07). (3) 
The distribution of FWHM([0 III] C ) reflects that NLS1 galaxies have similar [O III] C widths as 
BLS1 galaxies (Figure 2d). However, if we exclude galaxies with high velocity shifts (i.e., f[om] c 
> 75 km s _1 ), NLS1 and BLS1 galaxies show different scatter (the inset of Figure 2d) (KX07). 
(4) The distribution of velocity shifts of [O III] core lines (Figure 2e) reflects that NLS1 and BLS1 
galaxies differ in [O III] outflow velocity, in the sense that NLS1 galaxies have higher blueshifts of 
[O III] core lines than BLS1 galaxies, including a number of "blue outliers" with extreme blueshifts 
(K08). (5) The distributions of optical Fe II and [O III] strength (Figure 2f and 2g) indicate that 
NLS1 galaxies have stronger Fe II emission than BLS1 galaxies, and similar [O III] emission as 
BLS1 galaxies (X07). 

We list in Table 1 the key properties of our sample. Column (1) gives celestial coordinates in 
RA and DEC. Column (2) gives galaxy name. Column (3) gives the redshift. Column (4) lists the 
SDSS i-band absolute magnitude, derived from the PSF i magnitude. Column (5) lists the FWHM 
of the broad component of H/3 in km s _1 . Column (6) lists the ratio of total [O III] over total H/3 
emission. Column (7) lists the ratio of Fe II 4570 over total H/3 emission. Column (8) lists the 
FWHM of [S II] in km s^ 1 . Column (9) lists the FWHM of the core of [O III] in km s _1 . Column 
(10) lists the [O III] C velocity shift (blueshift) with respect to [S II] in km s _1 . Column (11) lists 
the [S II] intensity ratio of [S II] A6716/A6716. Column (12) lists the NLR electron density in cm" 3 , 
estimated from the density-sensitive [S II] intensity ratio (see X07 for details). Column (13) lists 
the flux density at 5100 A restframe in 10 -17 erg s _1 cm~ 2 A -1 . Column (14) lists the log of the 
monochromatic luminosity at 5100 A. Column (15) lists the log of the Eddington ratio. Column 
(16) lists the log of the black hole mass. The key sample properties of NLS1 and BLS1 galaxies are 
summarized in Table 2. 



4. Correlation analyses 

We now focus on studying trends across our sample, in order to perform a comparison with 
previously known correlations derived for different NLS1 samples, and to identify new trends. 
We have performed two types of correlation analyses. First, a Spearman rank order correlation 



analysis was used to determine which parameters are correlated and to derive the significance of 
the correlations. Second, we have performed a principle component analysis (PCA; e.g. Pearson 
1901, Duda et al. 2001). The PCA is a commonly applied method aimed at identifying the 
strongest correlations and the underlying parameters. This method allows us to link together diverse 
correlations into meaningful groups and reveal the physical mechanisms driving these correlations 
(for further discussion, and some caveats, see Boroson 2004). 

4.1. Spearman rank order correlations 

In a first step we computed Spearman rank order correlation coefficients between 11 of the 
key properties of our sample which were given in Table 1. These correlations involve emission-line 
widths (FWHM(H/3 b ), FWHM([S II]) and FWHM([0 III] C )), emission-line ratios (R5007, R4570, 
the ratio Rrsui defined as the intensity ratio of [S II] A6716/A6731, and the inferred NLR density 
n e ), and the parameters AL5100, L/LEdd, and Mbh- 

We note in passing that we have included both parameters, R[sn] and n e , in the correlation 
analysis. Even though strongly dependent on each other, their dependence is not linear, and we 
have therefore kept both parameters. Correlations with n c are slightly stronger than with Rpirj- 

Table 3 shows the results of our correlation analysis. We have split our sample into three 
groups; the NLS1 galaxies (55 objects), the BLS1 galaxies (39 objects), and the whole sample, 
and we have run the correlation analysis separately for each of the three, and separately report 
correlation coefficients in Table 3. All correlations with Spearman rank probabilities P s < 0.01 are 
marked in boldface. Figure 3 displays the correlation diagrams for several key parameters. 

The strongest (anti-)correlation among independent parameters involves the strength of Fe II 
over total H/3 emission, R4570. Among directly observed parameters, the clearest anti-correlation is 
between R4570 and the width of the broad component of H/3, FWHM(H/3b). The trend is significant 
for the whole sample, with a Spearman rank order correlation coefficient of r s = —0.7. This anti- 
correlation remains strong among the NLS1 and BLS1 galaxy sub-populations. In addition, we 
also confirm the known anti-correlation between the strengths of Fe II (R4570) and [O III] (R5007), 
which again is apparent in the full sample, and the two subsamples. The third strong correlation 
involving R4570, which persists across the BLS1 and NLS1 subsample, is the one between R4570 
and L/L^dd ( r s = 0.7). That correlation is stronger than with luminosity (no correlation at all) or 
with black hole mass, suggesting that L/LEdd is the more fundamental parameter. 

Two further correlations which involve the width of H/3b echo the trend we have already 
reported previously: the correlation between FWHM(H/3b) and density, which becomes apparent 
for the full sample (X07), and the correlation between W[oni] c and FWHM([0 III] C ) (K08), which 
is prominent in NLS1 galaxies (r s = 0.6), but absent in BLS1 galaxies. 

We have also investigated correlations among the widths of H/3b, [S II], and [O III]. While 



there is no correlation between FWHM(H/?b) and FWHM([S II]) for each subsample, such a trend 
becomes apparent across the full sample. The finding that a correlation between FWHM([S II]) 
and FWHMQO III] C ) is significant among BLS1 galaxies with r s = 0.67, but does not exist at all 
among the NLS1 galaxies, was already explained by K08: a subsample of NLS1 galaxies (so called 
blue-outliers) have a peculiarly broadened extra-component in [O III] (K08), and these few blue 
outliers destroy the expected correlation. If we remove from our NLS1 sample all those galaxies 
which show evidence for blue outliers (explicitly, objects with [O III] blueshifts exceeding i>[oni] c 
> 75 km s _1 ), then we also find a correlation between FWHM([S II]) and FWHMQO III] C ) in our 
NLS1 sample (r s = 0.47). 

We note that the width of H/?b does not systematically increase with continuum luminosity 
AL5100- This shows that NLS1 classification does not strongly depend on luminosity (Figure 4a). 
We also note in passing, that all of our NLS1 galaxies fulfill the criterion of Netzer & Trakhtenbrot 
(2007) of L/L Edd > 0.25 (Figure 4b). 

We would like to recall here the well-known fact that correlation analyses depend on the 
parameters considered, and they may not reveal the ultimate underlying correlation or driver. 
Two parameters may simply correlate because both depend on a third parameter which was not 
considered in the analysis. However, simple correlation analysis is a useful first step in searching for 
trends across data sets, and is a commonly used method, and in that sense we also applied it. As 
next step toward the final goal of identifying the strongest correlations in our sample and finding 
the key physical drivers behind them, we have applied a principal component analysis. 



PC A is a mathematical technique to convert a set of observations of possibly correlated vari- 
ables into a set of values of uncorrelated variables called principal components, or eigenvectors. 
The eigenvectors as a whole reproduce the original data space, and are orthogonal to each other. 
Eigenvector 1 (hereafter EV1) accounts for the largest variance of the data among the original mea- 
sured variables, and Eigenvector 2 (hereafter EV2), being orthogonal to (uncorrelated with) the 
first, represents the second greatest variance. Eigenvalues are defined as the fraction of the variance 
that each eigenvector accounts for. A good account of the application of PCA in astronomy can be 
found in Francis & Wills (1999). 

In a second step, we have run a pca| on the following input parameters: Mi, FWHM(H/?b), 
R5007, R4570, FWHM([S II]), v, om . c arid R[sn]- Each parameter represents independent informa- 
tion. Note that for the first time we involve the NLR density (reflected by the ratio Rfsui) hi the 
PCA. The results of the PCA are summarized in Table 4 (see also Figure 5 and Figure 6). 



4.2. Principal component analysis (PCA) 




We find that EV1 accounts for 38% of the variance in the data. EV1 is significantly (anti- 
)correlated with R4570, FWHM(H/3 b ), R [SII] , FWHM([S II]) and v [oul]c . When EV1 decreases, 
R4570 increases, outflow velocity becomes stronger and Rpii] increases, while R5007, FWHM(H/?b) 
and FWHM([S II]) decrease. EV2 accounts for 18% of the variance. EV2 is mostly anti-correlated 
with the absolute i magnitude and correlated with FWHMQS II]). 

In order to understand the relationship between the physical properties of AGN and the first 
two prime eigenvectors, we test correlations between EV1 and EV2 and n c , AL5100, Mbh, and 
L/LEdd, respectively (Figure 5), using the Spearman rank correlation test (Table 5). It can be seen 
that EV1 of our sample strongly correlates with n e , in addition to L/LEdd an d Mbh, while EV2 
with AL5100 and Mbh- 

Figure 6 shows the distribution of the objects of our sample in the EV1-EV2 space. In order to 
examine the dependence of the distribution on different parameters, each parameter is divided into 
three bins and coded by symbol size. The size of each point refers to coding according to [O III] C 
blueshift, Rjsii], FWHMjgn], R4570, L/Lgdd) AfBH, respectively. It is interesting to note that some 
of the objects with extreme values of EV1 and EV2 can be isolated from the EV1-EV2 diagram. 
We find that several "blue outliers" are located at the upper left corner of the figure (Figure 6a) . 
This regime is also dominated by objects with strong Fe II emission and high L/LEdd (Figure 6d 
and 6e). We also note that objects with high R[sn] fall at the low end of EV1 of the diagram 
(Figure 6b), and objects with high black hole mass are found preferentially at the high end of EV1 
(Figure 6f). 

5. The large-scale environment 

In order to probe the large-scale environment of the galaxies in our sample, and to search for 
possible differences in the environment of the NLS1 galaxies in comparison to the BLS1 galaxies, 
we have used the projected redshift-space two-point cross-correlation function (2PCCF), w p (r p ). 
We calculated the 2PCCF between the NLS1 or BLS1 galaxy samples and a reference sample of 
about half a million galaxies selected from the main spectroscopic sample of the SDSS final data 
release (DR7; Abazajian et al. 2009). A random sample was constructed so as to have the same 
selection effects as the reference sample. The reference and random samples were cross-correlated 
with the same set of galaxies (the NLS1 or BLS1 galaxy samples), and w p {r p ) as a function of 
the projected separation r p was defined by the ratio of the two pair counts minus one. Details 
about our methodology for computing the correlation functions and for constructing the reference 
and random samples can be found in Li et al. (2006). For consistency, we restricted our reference 
galaxies to the same redshift range as our NLS1/BLS1 samples, i.e., z < 0.3 

The amplitude of 2PCCF on scales larger than a few Mpc provides a direct measure of the mass 
of the dark matter halos that host the galaxies through the halo mass-bias relation. As shown in 
Li et al. (2008), the amplitude of the correlation function on scales < 100 kpc can serve as a probe 



of physical processes such as mergers or galaxy-galaxy interactions. On intermediate scales, the 
correlation probes the so-called '1-halo' term where the pair counts are dominated by galaxy pairs 
in the same dark matter halo. Therefore, the 2PPCF is a powerful measure of "environment" in 
the sense that it encapsulates information about how galaxy properties depend on the overdensity 
of the galaxy environment over a wide range of physical scales. 

In Figure 7 we compare the 2PCCFs estimated for our NLS1 and BLS1 galaxies. We do not see 
significant difference in w p (r p ) between the two classes of galaxies, over all the spatial scales probed 
(Wkpc/h < r p < 30Mpc/h). This indicates that, to the first order, NLS1 and BLS1 galaxies are 
found in similar environments PI 



6. Discussion 

Previous observations and their interpretation strongly indicate that NLS1 galaxies are AGN 
with low mass black holes typically accreting near the Eddington rate. As such, NLS1 galaxies as 
a class may hold important clues on black hole growth and evolution (e.g., Mathur et al. 2001), 
and on feeding and feedback in the course of galaxy evolution. We have therefore studied and 
homogeneously analyzed a large sample of NLS1 galaxies, accompanied by a comparison sample of 
BLS1 galaxies analyzed in the same way. This paper of the sequence has focused on correlation 
analyses, and PCA in particular. 



6.1. PCA and Eigenvector space 

PCA has previously turned out to be a powerful tool in uncovering trends and correlations 
among AGN samples (e.g., Boroson k Green 1992, Brandt 1999, Sulentic et al. 2000, 2002, Marziani 
et al. 2001, Grupe 2004, Wang et al. 2006, Sulentic et al. 2007, Zamfir et al. 2008, Mao et al. 
2009, Kuraszkiewicz et al. 2009, Marziani et al. 2010; see Boroson 2004 for a critical assessment 
of advantages and shortcomings of this method). However, it has been also noted that PCA are 
specific to the samples examined, depending on the observed properties used, as well as on the 
ranges of the parameters (e.g., Grupe 2004). Different eigenvectors may emerge from different 
samples. Adding new elements to PCA is therefore of great importance to verify the reality of 
the correlations, and furthermore to unveil the physical properties which determine the observable 
characteristics of AGN. 

The most popular interpretation of EV1 by far has been that EV1 is likely to be driven by 
the Eddington ratio L/L E dd (e.g., Sulentic et al. 2000, Boroson 2002, Grupe 2004). For the first 



3 It is well known that the clustering of galaxies strongly depends on their host properties. In order to perform a 
more detailed comparison, than we have done so far, a close matching of all involved samples in host galaxy properties, 
like luminosity, color and morphology, is required. 



time, we included the NLR density in the PCA, and have shown that EV1 of our sample is highly- 
correlated with density, in addition to L/Ledd- EV2 of our sample is dominated by luminosity and 
therefore accretion rate (proportional to luminosity), which matches the interpretation of EV2 of 
Boroson (2002). 

One key advantage of PCA is to discriminate between the various classes of objects according 
to their loci on the EV1-EV2 plane. The identification of objects that represent extremes are 
particularly crucial because their properties may suggest the nature of the physical parameter that 
governs the correlations (e.g., Boroson 2004). The similarities of the first two eigenvectors of our 
sample and Boroson (2002) allow a comparison of extreme objects in the EV1-EV2 diagram. 

We earlier found that NLS1 galaxies have higher blueshifts of the cores of their [O III] emission 
lines than BLS1 galaxies, including a number of "blue outliers" with extreme outflows (K08). 
Several of these blue outliers are located in the high Eddington ratio and high luminosity corner 
of the EV1-EV2 diagram (Figure 6). Interestingly, broad absorption line (BAL) QSOs, which 
show high outflow velocities in their absorption-line spectra (e.g., Weymann et al. 1991), fall into a 
similarly extreme locus (e.g., Boroson 2002) as the blue outliers of our sample, therefore pointing to 
possible links between these two source classes. Possibly, they bridge a gap between NLS1 galaxies 
and BAL QSOs (see the discussion of, e.g., Brandt & Gallagher 2000 on general similarities between 
NLS1 galaxies and BAL quasars). UV spectroscopy is needed to see whether the BAL fraction is 
high in these extreme blue outliers. 

6.2. Links between central engine and host properties 

Based on the correlation analyses we carried out, the strongest correlations among the directly 
measured parameters of our sample are between FWHM(H/3b), R4570 and R5007, confirming previ- 
ous studies (e.g., Boroson & Green 1992, Grupe 2004). This confirmation is of interest, because few 
independent large NLS1-BLS1 samples exist, for which a rigorous correlation analysis has been per- 
formed. Samples had different selection effects (e.g., soft X-ray selection), and therefore new NLS1 
samples which were selected differently from previous ones are important to check the persistency 
of trends. 

More importantly, we have shown that density is as important an ingredient in correlation 
space as Eddington ratio. Including density is of particular interest, because the NLR density is 
representative of the interstellar medium (ISM) of the host galaxy, and this study therefore links 
the NLS1 central engine properties with the host galaxy. Host properties are otherwise not easily 
accessible, and have only been studied for a few samples of NLS1 galaxies, focusing on the host 
morphology (Crenshaw et al. 2003, Deo et al. 2006, Ohta et al. 2007). Especially, host properties 
of NLS1 galaxies are an independent indicator of their evolutionary stage (Mathur et al. 2011, 
Orban de Xivry et al. 2011). 

Several factors can affect the density of the ISM of the host galaxy. Winds and outflows on 



the one hand (e.g., Schiano 1986, Kaiser et al. 2000), (bar-driven) inflows on the other hand (e.g., 
Shlosman et al. 1990, Riffel et al. 2008). We comment on each of them in turn. 

Winds/outflows are seen in nearby Seyfert galaxies (e.g., Kraemer et al. 2008, Wang et al. 
2010) and in NLS1 galaxies in particular (e.g., Komossa et al. 2008, Zhang et al. 2011). These can 
be disk-wind driven, or also be powered by radio jets (Morganti et al. 2010). The latter effect is 
unlikely to dominate in NLS1 galaxies. Even though a small fraction of them shows evidence for 
relativistic jets in radio and gamma-rays (e.g., Zhou et al. 2003, Yuan et al. 2008, Komossa et al. 
2006, Abdo et al. 2009, Foschini 2011), the majority of them is less radio-loud, on average, than 
BLS1 galaxies (Komossa et al. 2006), and rarely shows widely extended radio emission. Eddington- 
driven winds are more likely (e.g., Grupe 2004), given the high values of L/L^dd of NLS1 galaxies. 
These winds may extend into the inner NLR (Proga et al. 2008), and possibly further. Galaxy 
merger simulations predict large-scale outflows (e.g., di Matteo et al. 2005), but the merger fraction 
among NLS1 galaxies is still poorly constrained (e.g., Krongold et al. 2001), and we do not find a 
strong excess of pairing among our sample. 

On the other hand, recent studies indicate an excess of bars in NLS1 galaxies (Crenshaw et 
al. 2003, Ohta et al. 2007), and bar-driven instabilities might work at replenishing the NLR with 
(low-density) gas. Indeed, theoretical studies have shown that bars are efficient to transport large 
amounts of gas inward (e.g., Shlosman et al. 1990). Simulations have further shown that strong 
bars drive the formation of, and maintain grand-design nuclear dust spiral structures in the central 
kiloparsec of galaxies (e.g., Maciejewski 2004a,b). The gas loses its angular momentum to the 
bar via nuclear spiral shocks and is eventually concentrated in a narrow nuclear ring if it follows 
closely the periodic orbits in the bar (e.g., Patsis & Athanassoula 2000, Maciejewski et al. 2002). 
Regarding NLS1 galaxies, uninterrupted asymmetries able to drive gas inwards all the way from a 
few kpc down to a few tens of pc (e.g., Orban de Xivry et al. 2011) are detected. In particular, 
NLS1 galaxies show a higher fraction of grand-design dust spirals within ~ 1 kpc and stellar nuclear 
rings than BLS1 galaxies (Deo et al. 2006). Bar instabilities can also drive the formation of pseudo- 
bulges by internal secular processes (e.g. Kormendy & Kennicutt 2004). Such secular processes in 
NLS1 galaxy have been inferred based on bulge-disc decompositions (Ryan et al. 2007, Mathur et 
al. 2011, Orban de Xivry et al. 2011). 

Well-resolved host images are not yet available for the bulk of the NLS1 galaxies of our sample, 
for further distinguishing between the two scenarios discussed above, but could be obtained in the 
future with HST. Combining other host galaxy properties (structure, presence and properties of 
bars, stellar populations) with emission-line properties into correlation analyses will be another 
important future step. 



7. Summary and conclusions 



We have studied the correlations among the measured and derived properties of a sample of 
NLS1 and BLS1 galaxies on the basis of a principle component analysis. For the first time, we have 
involved the density of the narrow-line region (measured by the density-sensitive ratio of the two 
emission lines [S II] AA6716,6713) in such an analysis. A nearest-neighbor analysis of the large-scale 
environment of NLS1 galaxies was also performed. Our main results can be summarized as follows: 

• We have found that, among the parameters measured in this sample, the density of the 
narrow-line region is a key element of Eigenvector 1 (EV1), as important as the Eddington 
ratio L/LEdd- This is of particular interest, because it links the central engine and the host 
properties. 

• Apart from this new finding, we also confirm several previously known trends, especially the 
strong correlations involving the line width of H/3, and the strength of the Fe II complex and 
[O III] A5007 emission. In addition to density, L/L^d plays a significant role in affecting 
EV1 of our sample, while EV2 is related to luminosity. NLS1 and BLS1 galaxies are well 
separated in EV1 space, while they are not distinguished in EV2 space. 

• Accretion-driven winds on the one hand, or bar-driven inflows on the other hand, may play 
a role in explaining the links between the immediate vicinity of the supermassive black hole 
(SMBH) on the one hand (as traced by Eddington ratio and broad-line width, for instance), 
and the properties (especially density) of the host galaxy on the other hand. 

• Several galaxies with strong blueshifts of [O III] ("blue outliers") lie in the high Eddington 
ratio and high luminosity regime in the EV1-EV2 space, possibly sharing these properties with 
broad absorption line QSOs, and therefore suggesting possible connections between these two 
source classes. 

• To the first order, the NLS1 and BLS1 galaxies of our sample are found to reside in similar 
large-scale environments. While this needs to be confirmed with larger samples, it tentatively 
indicates the lack of an excess of mergers among NLS1 galaxies. 

Large-sample analyses of NLS1 galaxies and their remarkable properties are an important new 
approach in our understanding of black hole growth, accretion modes, feeding and feedback, and 
of aspects of galaxy - SMBH co-evolution and the role of secular evolution. 
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Table 1. Key properties of the NLS1 and BLS1 galaxies 
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025646.96+011349.4 


S0254+0101 





.177 


-22 


6 


2590 





.6 





.3 


570 


470 


30 


1 


.07 


450 


34.2 


44.0 


-0.5 


7.3 


030124.26+011023.0 


SDSSJ03014+0110 





072 


-20. 


7 


2940 





.4 


1 


.0 


250 


360 


110 


1 


.19 


250 


33.5 


43.3 


-0.8 


6.9 


030144.20+011530.9 


RXJ03017+0115 





075 


-21 


1 


3450 





.6 





.4 


340 


300 


40 


1 


13 


340 


43.2 


43.4 


-0.9 


7.2 


030417.78+002827.4 


KUG0301+002 





045 


-19 


7 


1200 





.6 





.1 


220 


110 


20 


1 


20 


240 


31.6 


42.9 


-0.2 


5.9 


030639.57+000343.2 


SDSSJ03066+0003 





107 


-21 


8 


1840 





.6 





.2 


330 


380 


-30 


1 


20 


240 


35.8 


43.6 


-0.3 


6.8 


031027.83-004950.8 


KUV03079-0101 





080 


-22 


4 


2820 





.3 





.5 


270 


130 


-10 


1 


19 


250 


150 


44.0 


-0.5 


7.4 


032213.90+005513.5 


KUV03197+0045 





185 


-23 


8 


2490 





.3 





.3 


300 


270 


50 


1 


.05 


490 


110.8 


44.6 


-0.2 


7.7 


032337.65+003555.7 


SDSSJ03236+0035 





215 


-22 


8 


1740 





.2 





.5 


190 


270 


110 


1 


.08 


430 


38.2 


44.3 


-0.04 


7.1 


032606.75+011429.9 


SDSSJ03261+0114 





127 


-21 


3 


1230 





.5 





.9 


330 


530 


180 


1 


.13 


340 


16.4 


43.5 


-0.01 


6.3 


032729.88-005958.5 


SDSSJ03274-0059 





134 


-21 


5 


5690 





.8 





.3 


290 


290 


20 


1 


.04 


510 


18.2 


43.5 


-1.3 


7.7 


034131.95-000933.0 


SDSSJ03415-0009 





223 


-21 


4 


1480 





.6 





.7 


200 


170 


10 


1 


.36 


60 


8.9 


43.7 


-0.1 


6.6 


075245.60+261735.8 


R.X .107527+ 26 17 





082 


-21 


4 


1600 





.3 





.8 


290 


210 


30 


1 


.06 


470 


55.1 


43.6 


-0.2 


6.6 


083949.65+484701.5 


NPM1G+48.0114 





039 


-19 


8 


1290 


1 


.2 





.2 


220 


190 


20 


1 


.19 


250 


50.4 


42.9 


-0.2 


6.0 


091313.73+365817.3 


RXJ09132+3658 





107 


-21 





1680 


1 


.0 





.5 


160 


350 


150 


1 


.38 


40 


21.3 


43.4 


-0.3 


6.5 


092247.03+512038.0 


SBS0919+515 





160 


-22 


3 


1250 





.3 


1 


.3 


160 


720 


430 


1 


.42 


10 


38.5 


44.0 


0.2 


6.7 


100405.00-003253.4 


SDSSJ10040-0032 





289 


-22 


6 


1410 





.4 


2 


.0 


150 


140 


40 


1 


.21 


230 


28.8 


44.4 


0.2 


7.0 


101645.11+421025.5 


RXJ10167+4210 





055 


-21 


2 


1690 





.2 





.8 


270 


250 


20 


1 


.23 


200 


92.4 


43.5 


-0.3 


6.6 


102434.72+555626.5 


SBS1021+561 





197 


-22 


4 


1660 





.6 





.7 


200 


240 


50 


1 


.20 


240 


29.9 


44.1 


-0.1 


7.0 


102448.57+003538.0 


SDSSJ10248+0035 





095 


-20 


8 


1990 





.6 





.3 


350 


230 





1 


.26 


160 


17.1 


43.2 


-0.5 


6.6 


102531.29+514034.9 


MARK142 





045 


-20 


6 


1530 





.2 





.9 


180 


250 


50 


1 


.28 


140 


97.2 


43.3 


-0.2 


6.4 


104210.03-001814.7 


SDSSJ10421-0018 





115 


-19 


9 


1150 





.3 





.9 


170 


210 


50 


1 


33 


90 


6.8 


43.0 


-0.1 


5.9 


104230.14+010223.7 


SDSSJ10425+0102 





115 


-20 


8 


1430 


1 


1 





.5 


260 


300 


90 


1 


.18 


270 


13.4 


43.3 


-0.2 


6.3 


104332.88+010108.8 


SDSSJ10435+0101 





072 


-20 


4 


2530 


1 


.0 





.1 


230 


200 





1 


.19 


250 


18.4 


43.0 


-0.8 


6.6 


112941.94+512050.7 


SBS1126+516 





234 


-23 


2 


1880 





.6 





.3 


190 


250 


60 


1 


.17 


280 


52.1 


44.5 


-0.04 


7.3 


115023.59+000839.1 


Q1147+0025 





127 


-21 


1 


1310 





.5 





.6 


210 


150 








.97 


690 


16.7 


43.5 


-0.1 


6.4 


115235.00-000542.8 


Q1150+0010 





129 


-21 


6 


3140 





.4 





.4 


280 


190 


10 


1 


.14 


330 


27.1 


43.7 


-0.7 


7.3 


115533 50+010730.6 


SDSSJ11555+0107 





197 


-21 


9 


1510 





.3 





.7 


150 


780 


330 


1 


.23 


200 


17.1 


43.8 


-0.05 


6.7 


115758.73-002220.8 


QUESTJ1157-0022 





260 


-23 


8 


4670 





.6 





.3 


360 


280 


-10 


1 


.20 


240 


93.4 


44.8 


-0.7 


8.4 


120226.76-012915.3 


IRAS11598-0112 





.150 


-22 


4 


1460 





.5 


2 


.7 


250 


340 


170 


1 


.31 


110 


47.6 


44.1 


0.1 


6.8 



Table 1 — Continued 



coordinates (J2000) 


common name 


z 


Mj 


FWHM(H/3 b ) 


R.5007 


R4570 


FWHM([S II]) 


FWHM([0 III] C J 


"[OIII]c 


R [SII] 


n e 


/5100 


log AL 51 nn 


log L/L Edd 


logM BH /M 




(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


121549.45 + 544224.0 


b-tib lzld + o4yA 


0. 150 


-22.9 


1240 


0.6 


1.2 


310 


320 


110 


1 


.27 


150 


73.4 


44.2 


0.3 


6.8 




OROCmO/l'J/l 1 A 1 A A 

blJbb.J lz4d4-|-UlUU 


0.090 


-20.6 


5180 


0.4 


0.2 


230 


otn 


-10 





.94 


IIS 


16.5 


Ti 


-1.3 


7.3 


1 O^KO'J AA I aaocja (1 
124ozd.UU-r-UUzfSdy.y 


OnOC 11 O/IC'J 1 AAOQ 
dUDDJ lz4Dd-|-UUZS 


0.088 


-20.6 


2 onn 


0.9 


0.3 






40 


1 


.14 


330 


17.6 




-0.6 


6.6 


i O/iRor, or, i aoooaq c 
1 Z4bdo . zo-(-UzzzUo. o 


Dfn/i/i 1 AO*; 
r\j 1 Z44-|-Uzo 


0.048 


-20.5 


1200 


0.6 


0.8 


230 


300 


150 


1 


.10 


400 


82.0 


43.3 


-0.03 


(i.2 


1 9AA19 no AAE/IOA « 

loUUzo.zz — UU04zy.o 


T T A H E •} A 

U lvlod4 




"00 t 


oo/tn 


0.7 


0.6 


250 


150 


70 


1 


.09 


ton 


20.9 


43.5 


0.01 


6.3 


1 OI175 f O E A A 'J CI A 1 « 

loU/ Id. ZD — UUooUl.o 


dUddJ loU I Z-UUoO 


0.170 


-22.5 


2240 


0.9 


0.3 


370 


310 


-30 


1 


.01 


■:>80 


40.6 


44.1 


-0.3 


7.2 


i »)T7rn O^l 1 A 1 A A 1 A 

Idl f DU.dz + o01041.0 


O ~D O 1 O 1 E 1 i?A/1 

b -tJ b 1 d 1 b + 04 


0. 137 


-22.1 


2660 


1.5 


0.1 


380 


340 


20 


1 


.20 


240 


31.3 


43.8 


-0.6 


7.2 


1 d Z 1 d O . d d — U U 1 d U . 


CnCC 11 ')OI r, A A 1 'J 
bUbb.J IdZlO-UUld 


0.082 




3300 


1.0 


0.3 


^t? 


300 


130 


1 


.11 


380 


38.5 


43.5 


-0.9 


7. 1 


1 O/IQCi (|7 1 AAA/1 'J /I C 

ld4d Jl .U I +UUU4d4.» 


CROC 11 Q/I'JB I A A A A 

bUbb.J ld4do-|-UUU4 


0.074 


"on 7 
"00 I 


1960 


0.7 


0.4 


200 


190 


10 





.99 


630 


2 ? t 


/M \ 


-0.4 


6.5 


1 O/IHCn /IE AA1 CCO E 

ld44&y.4 [ ) — UOlJJM.J 


<"l 1 'J. A O AAA 

l»Jld4z-UUU 


0.245 


-23.0 


2450 


1.0 


0.2 




Si/in 


-30 


1 


.26 


160 


41.7 




-0.3 


7.5 


loecii; EE i ci;i oi/l "7 
Id J JlO.O J + ODlz44. f 


b-Dbld jd + OD4 




"on 4 


1340 


1.6 


0.7 


490 




50 





.94 


770 


60.2 


aa n 


0.1 


6.7 


1 /lOAQA 11 AA1 1 1 E 1 

14dUdU . zz — UU 1 1 1 . 1 


CnCC T1 /lonr, A A 1 1 

bUbb.J 14dUo-UUl 1 


0. 103 




1510 


1.3 


0.5 


170 


160 


30 


1 


.23 


200 


8.4 


43.0 


-0.3 


6.2 


1 /I Q A O A "} 1 1 /lA'JEOA A 

14dUdy .d 1 -|-4ydody.U 


L/b*-JDDl 


0.204 


"00 7 




0.3 


0.5 


260 


230 


30 


1 


19 


2 sn 


38.9 


tt 


-0.4 


7.4 


1 AOQAT E A AAAQAE A 

14do4 i . o4 — UUUoUo .4 


CRCC 11 /l')U7 A A A Q 

bUbb.J 14do f-UUUo 


0.104 


"„„„ 
-20.9 


3570 


0.2 


0.4 


200 




-20 


1 


.01 








-1.0 


7.1 


1 1 1 1 OA 11 1 EAOQA1 7 

144idU. 1 1 -(-oyzoUl . 1 


DV111111 1 rflOO 

UAJ 14414-j-oyzo 


0. 134 




1600 


0.3 


0.8 


180 


250 


30 


1 


.36 


00 


16.8 


43.5 


-0.2 


6.6 


144 / oj.zd — UUdZdU . <-} 


y nc O T 1 A A 7e. AA'iO 

bUbb.J 144 ( O-UUdz 


o'osl 


21 6 














1 


.20 


?on 


Js'n 








1 4 4 ; } .. i Z . 1 1~|-UUZZOO.O 


bubb.i i44yo-r~uuzz 




"on 5 


104f! 


7 


n 8 


230 


lin 


in 


1 


.30 






43 1 


1 


6 1 


±401Zd.UZ — U l J l J t ) Z i ) . : J 


encq ti a r.i 'J nnnfi 

bUbb.I 14Old-UUU0 


1 30 


21 1 




8 


1 


200 




" 


1 


.23 


200 




43 4 






1 /I K1 /) Q OA 1 KO/11 17 /I 
14D 14d . dU~rOZ41Z ( .4 






"00 n 






no 








1 


.04 




18 3 




"(13 




1 1 CCQi i5E AA1 1 1 A O 

14oDol.DO — UU 1 1 14. Z 


OnQO T1 ircr A A 1 1 
blJbb.J 14O0O-UU 1 1 


n 1 


"01 t 


roon 


0.5 


0. 1 




o,n 

230 


on 





.87 


1030 


25.3 




" , 
-1.3 


7.8 


1 Ei QEiS C7 I (1111 C1 /I « 


CPiot; T1 ^1 oo_i_AAif; 


114 




7-)n 






n 


?7n 


"an 


1 


.07 


fso 




43 6 






i f,of;os on nmssnQ a 

lOZOZo.zU — UUdoUy .4 


bUbb.I loZ04-UUdo 


123 


21 5 




3 


n r 






20 


1 


.25 




27 4 




"n 5 


6 Q 


i E0701 fiO-L/l Q/1 9/1 7 s 
1 Jd / dz .DZ-j-4y4Z4 ( . 


QRQl ^"if^-L/l OSi 

b.t>blJdO-|~4yo 


280 


22 2 




5 


10 


120 


14(1 


30 


1 


.29 




21 6 


44 2 


"0 2 




lOijyii.i / -j-uuzouu. « 


ervee t i con 1 1 nnoa 

bUbb.I lOdyi-T"UUZ0 










4 






1 in 


1 


.25 


180 


28 2 


44 3 


01 


7 n 


1EEQ00 1 A 1 Q7AQ Q Q A 

looyzz.iy+z (Uddo.y 


0TT1 f,r,7 1 171 
ZHilDDf -|-Z / 1Z 




"on 7 




5 


n't 








1 


.23 


200 




43 2 


(17 


6 7 


ioiouy.do-j-doiyo * . o 


itAj ioioi-j-doiy 




1Q 1 






n q 






"on 


1 


.43 


< 10 






"no 




IfilAEI Ol 1 /1AEO/17 

lt)iyol.dl +4Uuo4 ( . z 


Ti'TT/^lfilQ 1 /I1A 
lS.UtjlOlo-|-41U 


0.038 


1Q O 

0? 7 


7cn 

17o0 


1 


0.2 


280 


on 


on 


1 


.16 


300 


oo'r 

28.5 


42.7 


" ,V 
-0.6 


6.1 


1 ^1901 CQ 1 /I71 OCO A 
lt)ddzd.OO+4 1 lOOy.U 


DYIIfilll 1 4 71 a 

r\.A.J lDddd+4 ( lo 


0.116 


-21.7 


1260 



1.2 


0.9 


340 


250 


An 
40 


1 


.17 


280 


38.0 


43.7 


0.1 


6.5 


lt)4yU / .04+o4z4zZ.d 


t~\ 00 n Cyim i ^/io/i 
bJJbbJ lD4yi + t)424 


0. 184 


-21.1 


1360 


0.5 


0.8 


160 


320 


10 


1 


.41 


20 


9.4 


43.5 


-0.1 


6.4 


165408.16+392533.4 


EX01652. 4+3930 


0.069 


-20.4 


1280 


0.2 


0.6 


250 


120 


50 


1 


.11 


380 


25.1 


43.1 


-0.2 


6.1 


165658.38+630051.1 


SDSSJ16569+6300 


0.169 


-21.1 


1720 


0.2 


0.7 


160 


210 


20 


1 


.15 


310 


10.8 


43.5 


-0.3 


6.6 


170328.97+614110.0 


SDSSJ17034+6141 


0.077 


-21.5 


4000 


0.6 


0.2 


260 


430 


60 


1 


.16 


300 


55.1 


43.6 


-1.0 


7.4 


170546.91+631059.1 


SDSSJ17057+6310 


0.119 


-20.9 


1930 


0.2 


0.7 


250 


230 


80 


1 


.37 


50 


14.9 


43.4 


-0.4 


6.6 


170812.29+601512.6 


SDSSJ17082+6015 


0.145 


-21.0 


1070 


1.0 


0.5 


280 


260 


60 


1 


.25 


180 


10.8 


43.4 


0.1 


6.1 


171411.63+575834.1 


SBS1713+580 


0.093 


-22.2 


2130 


0.4 


0.7 


360 


300 


6 


1 


.13 


340 


88.7 


43.9 


-0.3 


7.1 


171550.49+593548.8 


SDSSJ17158+5935 


0.066 


-20.0 


3240 


0.7 


0.2 


340 


270 


40 


1 


.14 


330 


15.8 


42.9 


-1.0 


6.7 


171829.01+573422.4 


SDSSJ17184+5734 


0.101 


-21.3 


1760 


0.4 


0.7 


110 


470 


150 


1 


.40 


30 


26.9 


43.5 


-0.3 


6.6 


172032.29+551330.3 


SDSSJ17205+5513 


0.273 


-22.7 


3490 


0.6 


0.2 


250 


250 


-10 


1 


.20 


240 


33.7 


44.4 


-0.6 


7.8 


172533.07+571645.6 


SDSSJ17255+5716 


0.066 


-20.5 


4870 


1.4 


0.04 


280 


210 


20 


1 


.26 


160 


23.1 


43.0 


-1.3 


7.2 


173107.87+620026.1 


SDSSJ17311+6200 


0.069 


-20.8 


3860 


1.1 


0.4 


350 


310 


20 


1 


.07 


450 


30.5 


43.2 


-1.1 


7.1 


221918.53+120753.2 


IIZwl77 


0.081 


-21.3 


1180 


0.5 


1.2 


170 


200 


40 


1 


.29 


130 


46.6 


43.5 


0.1 


6.3 


232328.00+002032.9 


SDSSJ23234+0020 


0.120 


-21.2 


3890 


0.5 


0.3 


410 


300 


10 


1 


.11 


380 


18.4 


43.5 


-1.0 


7.3 


233032.95+000026.4 


SDSSJ23305+0000 


0.123 


-21.3 


1220 


0.7 


0.8 


410 


280 


80 


1 


.15 


310 


18.5 


43.5 


0.01 


6.3 


234725.30-010643.7 


SDSSJ23474-0106 


0.182 


-21.8 


1830 


0.3 


0.4 


320 


180 


40 


1 


21 


230 


18.2 


43.8 


-0.2 


6.9 



a columns from left to right: (1) coordinates in RA (h,m,s) and DEC (d,m,s), (2) galaxy name, (3) redshift, (4) SDSS absolute i magnitude, (5) FWHM of the broad component of H/3 in km s 1 , (6) ratio 
of total [O III] over total H/3 emission, (7) ratio of Fell 4570 over total H/3 emission. (8) FWHM of [S II] in km s -1 , (9) FWHM of the core of [O III] in km s _1 , (10) velocity shift (blueshift) of the core of 
[O III] with respect to [S II] in km s _1 . (11) intensity ratio of [S II] A6716/A6716, (12) NLR electron density in cm" 3 , (13) flux density at 5100A restframc in 10~ 17 erg s -1 cm -2 A -1 , (14) log of the 
monochromatic luminosity at 5100 A, (15) log of the Eddington ratio, (16) log of the black hole mass. 



Table 2. Sample properties 



Prop. 


class 


range 


average 


log Mbu/Mq 


NLS1 


[5.7,7.3] 


6.5 




BLS1 


[6.5,8.4] 


7.2 


log L/L Edd 


NLS1 


[-0.6,0.3] 


-0.1 




BLS1 


[-1.4,-0.3] 


-0.8 


FWHM([SII]) 


NLS1 


[110,490] 


230 




BLS1 


[120,570] 


300 


FWHM([OIII] c ) 


NLS1 


[110,780] 


280 




NLSl a 


[110,440] 


210 




BLS1 


[140,470] 


270 




BLSl a 


[140,470] 


270 


R4570 


NLS1 


[0.1,2.7] 


0.7 




BLS1 


[0.0,1.0] 


0.3 


R5007 


NLS1 


[0.2,1.6] 


0.6 




BLS1 


[0.2,1.8] 


0.7 


U [OIII]c 


NLS1 


[-30,430] 


80 




BLS1 


[-60,130] 


10 


R[SII] 


NLS1 


[0.94,1.43] 


1.22 




BLS1 


[0.87,1.47] 


1.13 


n e 


NLS1 


[10,770] b 


200 




BLS1 


[150,1030] 


340 



a Excluding galaxies with high velocity shifts 

(i- e -' Vn]c > 75 km s_1 )- 

b Two objects, RXJ16181+3619 and 
SDSSJ01148-0029, have n e < 10cm" 3 , and 
were not included in the analysis involving 
n e . 



Table 3. Spearman rank order correlation a 



Prop. 


FWHM H/3k 


R5007 


R4570 




FWHMroilll 

[Wllljc 


V 

[OIIIJc 


Rrsui 

[Oil] 


n c 


AX/5100 


L/ ijE dd 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 
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a All correlations with Spearman rank probabilities P s < 0.01 are marked in boldface. 



Table 4. Results of the PCA. The relative significances of the eigenvectors are listed, as well as 

the projections of the original parameters. 
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1.2295 
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0.6432 
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0.4366 


0.3966 


Percentage variance 
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9.19 
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6.24 


5.67 


Cumulative 


37.99 


55.55 
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88.09 


94.33 


100.00 
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-0.885 


0.183 


-0.134 


-0.295 


0.259 


-0.020 
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0.277 


0.390 


FWHM [SII] 
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-0.304 
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0.025 
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R[sn] 
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Table 5. Spearman rank order correlation for the first two eigenvectors with n e , AL5100, Mbu 

and L/L Edd a 





EV1 


EV2 


n e 


0.71 


0.05 


AislOO 


-0.01 


0.86 


M BH 


0.60 


0.48 


L/ -^Edd 


-0.71 


0.36 



a All correlations with 
Spearman rank probabili- 
ties P s < 0.01 are marked 
in boldface. 
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Fig. 1. — Distribution of black hole masses and Eddington ratios of the galaxies of our sample. 
NLS1 galaxies are shown as solid line, BLS1 galaxies as dashed line. 
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Fig. 2. — Histograms showing the distribution of the sample properties of NLS1 and BLS1 galaxies. 
NLS1 properties are shown as solid line, BLS1 properties as dashed line. Panels from left to right 
and top to bottom: (a) FWHMQS II]) in km s -1 ; (b) intensity ratio of [SII]A6716/A6731; (c) 
logarithm of NLR electron density in cm -3 . The shaded histogram indicates the two galaxies with 
n e < 10cm -3 . (d) FWHMQO III] c ) in km s _1 . The inset shows the distribution of [O III] C velocity 
(blueshift) after excluding objects in the highest velocity bin; (e) [O III] C velocity (blueshift) w.r.t. 
[S II] in kms -1 ; (f) ratio of Fe II 4570 over total H/3 emission; (g) ratio of total [O III] over total 
H/3 emission. 
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Fig. 3. — Strongest correlations among emission- line and continuum properties. NLS1 galaxies are 
coded as filled circles and BLS1 galaxies as open circles. The objects that are off the plots are 
indicated by arrows. 
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Fig. 4. — FWHM of the broad component of H/3 versus AL5100 and L/LEdd- Symbols are as in 
Fig. 3. The dashed line in the right panel corresponds to L/L^dd = 0.25. 
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Fig. 5. — Correlations of the first two eigenvectors with L/LEdd> -^bh, n e , and ALsioo- Symbols 
are as in Fig. 3. The eigenvectors are plotted in arbitrary units. 
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Fig. 6. — Influence of different parameters on the distribution of NLS1 galaxies (filled circles) and 
BLS1 galaxies (open circles) with respect to the first two eigenvectors. Each parameter is divided 
into three bins and coded by circle size: large circles, high value; medium circles, intermediate value; 
small circles, low value. Left to right and top to bottom: coding according to [O III] C blueshift (blue 
outliers in [O III] C are marked with an extra open square.), Rpgn], FWHM([SII]), R4570, L/L-^^, 
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Fig. 7. — Projected cross-correlation function w p (r p ) for NLS1 galaxies (filled circles connected by 
the solid line) and BLS1 galaxies (open circles connected by the dashed line) in our samples, with 
respect to a reference sample of about half a million galaxies in the SDSS/DR7. Errors in the 
w p (r p ) measurements are estimated from the Bootstrap resampling technique. 



